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The fractional Fokker-Planck equation, which contains a variable diffusion coefficient, is discussed 
and solved. It corresponds to the Levy flights in a nonhomogeneous medium. For the case with the 
linear drift, the solution is stationary in the long-time limit and it represents the Levy process with 
a simple scaling. The solution for the drift term in the form Asgn(a;) possesses two different scales 
which correspond to the Levy indexes fi and /i + 1 (fj, < 1). The former component of the solution 
prevails at large distances but it diminishes with time for a given x. The fractional moments, as a 
function of time, are calculated. They rise with time and the rate of this growth increases with A. 
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Transport processes in physical systems can be very complex. The traditional statistical description, which relies 
on uniform Gaussian noises and position-independent transition probabilities, must fail for many realistic problems. 
A need for an alternative approach is obvious, for example, when one considers porous, irregular materials, containing 
impurities or entanglements which act as obstacles and dynamical traps Transport in a heterogeneous, in particular 
fractal, material involves position-dependent transition rates and highly complex driving forces which can be handled 
in a statistical manner (quenched disordered media) Q. The Gaussian distribution not always applies. In systems 
which are characterized by long-range correlations and non-local interactions, one can expect the presence of long tails 
of the driving noise, i.e. one should consider the general Levy distributions. It is so for many physical phenomena 
[H , in particular in biological 0, Q , social @ , epidemiological problems 0] ■ From the Langevin equation, driven by 
the homogeneous Levy noise, follows the Fokker-Planck equation (FPE), which is fractional @,@) E3, Ell- The 
importance of the general Levy distribution stems from its stability: it acts as an attractor in the functional space 
and there are no other attractors. The physical reason behind the Levy, non-Gaussian, processes traces back to 
the nonhomogeneous structure of the environment, in particular fractal or multifractal. However, this basis and the 
essential feature of the Levy process is only rarely taken explicitly into account and transport processes are described 
in terms of linear stochastic equations. The medium structure is usually reflected only in the form of the external 
potential or as a time-dependence of the stochastic driving. In this paper, we consider a nonhomogeneous Levy noise 
which leads to a position-dependent diffusion coefficient in FPE and construct the asymptotic solution of FPE. 

The problem of nonhomogeneous stochastic driving can be posed in a form of the Langevin equation, 

dx(t) = F{x)dt + a{x)dL, (1) 

with the multiplicative noise which is understood in a sense of the Ito interpretation. Eq.JT]) can be regarded as a 
result of the adiabatic elimination of fast variables in nonlinear processes with additive fluctuations. In the Gaussian 
case, it is suited e.g. for the following problems: the ensemble of two-level atoms in the electromagnetic field (Maxwell- 
Bloch equations), the parametric generation of coherent fields by incoming laser field, the Raman scattering, and the 
autocatalytic reactions [13j |. In the present paper, we assume that the noise L(t) is the uncorrelated Levy process 
with the stability index fi (0 < fi < 2), the median 7, and the scale parameter K^. The cumulant expansion of the 



characteristic function, truncated at the order jj,, produces the following fractional FPE [14 1 



^p(x,t) = - A [fra( x )+F(x))p(x,t)]+K»-^[D{x)p(x,t)] (2) 

for the probability density distribution of the variable x; in the above equation d^/d\x\^ is the Riesz fractional 
derivative and D(x) — |a(a;)|'*. The initial condition is p(x,0) = S(x). Eq.([2]) has been extensively studied for 
D(x) =const fl5j |. 

On the other hand, equations in the form ([2|) follow directly from master equation for a jumping process: 

p(x,t) = / dx'[w(x\x')p(x' , t) — w(x'\x)p(x, t)]. (3) 



For example, modelling the thermal activation of particles within the folded polymers leads to the FPE with the 
variable diffusion coefficient, which results from the polymer heterogeneity [16j |. In general, D{x) in FPE is variable 
if the transition probability w{x\x') ^ w{\x — x'\). It is the case also for a coupled continuous time random walk 



2 



(CTRW) which is defined in terms of a Poissonian waiting time distribution with a variable jumping rate v(x), as 
well as of the Levy distribution of jumping size Q(x) 17|, [l8|. Then w{x\x') = v(x')Q(\x — x'|); FPE takes the form 



(|2"|) with D(x) — v(x). The parameter 7 can depend, in general, on the position before the jump. If it is constant, the 
drift term takes the form ~ ypjx), where 7 = {x)q. The solution of that equation for the power-law v{x) and 7 = 
is still the Levy distribution [181 ]. 

In this paper we demonstrate that the stability property holds also for some systems which contain the driving 
term F(x). That is by no means obvious; for example, the separable solution of the fractional Schrodinger equation in 
Ref. [Ill , which is characterized by the exponential diffusion coefficient and a periodic potential, looses its dependence 
on [i altogether, in the asymptotic limit. For a system with a power-law external potential, |x| c , driven by the Levy 
noise with the constant diffusion coefficient [l9| , the stochastic properties depend on c: the Levy index in the stationary 
solution remains unchanged for the harmonic potential, whereas for larger powers the asymptotics is determined by 
c. As a result, the variance can be finite. We will solve Eq.© for the power-law diffusion coefficient, D(x) = \x\~ e 
[20I ] . and for two simplest forms of F(x) which correspond to symmetric potentials. 

It is convenient to handle Levy processes by means of the Fox functions. If the solution of the FPE is to be the 
Levy process, one can expect it has the following scaling form 



p(x,t) = Na(t)H$ 



a(i)|x| 



(ax,Ax), (a 2 ,A 2 ) 
{b 1 ,B l ),{b 2 ,B 2 ) 



(4) 



where N is the normalization constant. The expression © is universal for any Levy process 2l|. We assume the 
solution of Eq.© in this form. To find the coefficients, we require that Eq.© is satisfied by © in the diffusion (fluid) 
limit of small wave numbers k. Therefore, the solution in the form ((4]) should coincide with the exact solution of 
Eq.© for large \x\. 

We begin with the case F(x) — —Ax (A > 0) which corresponds to the harmonic oscillator potential. Moreover, we 
assume 7 = and fi + > 0. The FPE takes the form: 



d , . d , 

— p{x,t) = X—[xp(x,t)\ 

and the Fourier transformation produces the result 

— P (k,t) = -Xk—p{k,t) - 



<9 M 
d\xj 



9 p(x,t)} 



9 p(x,t)}. 



(•5) 



(6) 



To solve Eq.©, we follow the procedure from Refs. HI, 12], where also the appropriate formulas are provided. Due to 
the multiplication rule, the argument of the Fourier transform in the last term is the Fox function of the same order as 
p(x,t). We apply the formula for the cosine Fourier transform and expand the results in the fractional powers of \k\. 
In order to adjust the terms on both sides of Eq.©, we have to introduce conditions on the powers and to eliminate 
some of the terms by an appropriate choice of the coefficients. As a result, we determine the following coefficients of 
the Fox function: a x = 1 - (1 - 0)/(/x + 0), A 1 = 1/(^ + 0), b 2 = 1 - (1 - 0)/(2 + 0), and B 2 = 1/(2 + 0). Those values 
are - for any choice of the other parameters - a sufficient condition for Eq. (TJJ to represent the Levy process in the 
lowest order of the ^-expansion: 



p(k,t) ^il-hf, 



(7) 



For the Fourier transform in the last term in Eq.©, we need to keep only the term k . We obtain p~g = 
!F c [\x\~ p(x,t)] w h o a , where a results from the transformation to the scaled variable. The neglected terms on 
both sides of Eq.© are of the order |A;| 2ai+ . The expansion of the Fox function around zero and 00 shows that 
(61, Bi) corresponds to the behaviour of p(x,t) at x = 0, whereas (a\,Ai) determine the asymptotics (|x| — > 00). 
Therefore, the former ones cannot be determined in the small k approximation. We assume values which correspond 
to the small |x| limit of the master equation solution for CTRW [22|: b\ = and B\ = \; that process is described by 
Eq.© with A = 0. To settle a 2 and A 2l which only weakly influence p(x,t) in the asymptotic limit, we require that 
the x-dependence of the distribution © for A = should coincide with the stretched-Gaussian exact solution in the 



limit fj, -> 2 [22j; then: a 2 = 1/2+ (1 - 0)/(2 + 
the expansion formula: — N(p, + 0)T(— n)T(l 
whereas ho = J Q pg(ax,t)d(ax) — lim s _,_i^(s) 



and A 2 = (1 + 0)/(2 + 6 
-fi + 9) cos(7r^/2)/r[l/2- 
= AT(//+0)/(2 + 0)r[l/2- (0 2 



. The coefficient hp follows directly from 
(H + /i0 + 2)/(2 + 0)]r[-(/x + 0)/(2 + 0)], 
1 - 2)/(2 + 0)], where x(s) is the 
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Mellin transform from the Fox function. Similarly we obtain the normalization constant: N = r[— 0/(2 
2/(2 + 0)]/2r(l + 0)r[-0/(/x + 0)]. Inserting p and p e into Eq.© yields the equation for a(t): 



0)]r[i/2 



d = Xa - K* 1 



hp 



-a 



fj,+e+i 



which can be solved by separation of the variables: 



a(t) = 



X/c L 



1 - exp[-A(^ + 6)t] 



(8) 



(9) 



where cl = K^ho/ fih^. For = 0, the above solution agrees with that of Ref . (13] . 

Expansion of the Fox function in powers of l/\x\ reveals the asymptotics of the Levy process: p(x, t) ~ a(i)~ M |a;|~ AI ^ 1 
for |x| — > oo. a(t) approaches with time a constant which corresponds to the stationary solution of FPE. The speed 
of that convergence depends on 0: it is rapid for large 0, whereas negative values of can substantially hamper the 
convergence. The meaning of the parameter in the context of the diffusion process becomes clear when we consider 
the case for which the variance exists, namely the Gaussian case (x = 2; it is involved in solution ((4]). The coefficients 
of the Fox function on its main diagonal are then equal and wc can apply the reduction formula. The solution reads 



1,0 



p(x,t) = Na(t)H?" 



a(t)\x\ 



(i + 1 

V2 ^ 2 



1+0 ' 
2+9 ' 2+6' 



(M) 



(10) 



The contribution to the Barnes-Mellin integral from the residues vanish for large \x\ and the asymptotic form of the 
Fox function is stretched-exponential 13, 23l |: 



p[x,t) - 

where a(t) is given by Eq.© with fi = 2 and C2 



,1+6 



(1 + 0) 



exp[- C2 (a|x|) 2+e ], (11) 
2+6 1 (2 + 9) 3+e . The variance can be easily evaluated: 



x ^ = ~g^2 p { '^ = h2a 



(12) 



where \i2 = lim A1 ^2 hp. If A = 0, Eq. (fT2| predicts the normal diffusion (0 = 0), subdiffusion (0 > 0), and superdiffusion 
(0 < 0). Otherwise, the variance converges with time to a constant. Therefore, the parameter governs the transport 
speed. In the coupled CTRW, a large means that the average trapping time strongly rises with the distance. 

The system, which has been discussed above, is characterized by the same stability property as that of the driving 
noise: it is Levy distributed with the parameter /i. Moreover, it reveals the simple scaling. One can ask whether 
the same properties hold for other systems driven by the Levy noise and a symmetric potential 2J|. The next case 
demonstrates that, even if the stability property is preserved, the index n may change. Let us consider the drift in 
the form F(x) = Asgn(x), which corresponds to the wedge-shaped potential. We assume < p, < 1, i.e. a process of 
the infinite mean, and /i + > 1 . The FPE is the following 



d_ 

dt 



p(x,t) 



d 

A— [sgn(a;)p(a;,f)j 



d\xj 



9 p(x,t)]. 



Its cosine Fourier transform reads 
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dt 



p(k,t) 



- l p{x,t))- K»\k\»F c [\x\- e p{x,t)\ 



(13) 



(14) 



and the factor \x\ 1 , which results from the change of the sine to cosine transform, introduces a new scale. We take 
into account that double scaling by assuming the solution in the form p(x, t) — N[pi(x, t) 4- ap2(x 1 1)], where 



Pi{x,i) = fi{t)H. 



2,2 



.MOM 



(a®,A<p),(aS J ,W) 



(I, 



(15) 



and a is determined by the initial condition p(x, 0) = 6(x) — [S(x)+aS(x)]/ (1+a); we assume that a ^ if A ^ 0. The 
solving method is similar to the previous case: we insert Eq. (|15[) into Eq. (|14p and compare the terms. The essential 
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FIG. 1: (Colour online). The ratio which determines the relative contribution to the solution of Eg, (|13[l from the terms which 
correspond to the Levy indexes ji and [s, + 1, for various parameters of the process. All curves fall like 1/t. 



point is to realize that the first term on rhs upgrades the index \x by one, because the expansions of T c [\x\ 1 pi(x,t)) 
are determined by the terms |fc|( 2-a i ' . Then J 7 c [\x\~ 1 pi(x, t)] ss const — h'^ f^^\k\^ +1 : where terms of the order 

H + 2 have been neglected. We put — a,i, — A4, b^ — bi, and B 1 - 1 ' = Bi. The coefficients for p 2 are the same 
except fi — ► /j, + 1. By comparing the terms of order /j, and [i + 1, we obtain a set of two differential equations 



l _ K^r r- e/>i 1 n <-- 6/(fl+lh 



\h>v 



(/* + !) — -ferti, (16) 



ah 



where = /f, £ 2 = hfi = h„ = h^+x, = ho, h ( 2) = h (» -> M + 1), and /i'^ = -7V( M + 0)r(-// - 

l)r(2 + fi) sin(7r^/2)/r[l/2 + (ji + fj.6 + 2)/(2 + 6)]T[-([i + 0)/(2 + 9)}. 

Let us assume A > 0. The asymptotic form of p(x, t) involves contributions from both fi and fj, + 1, p(x,t) = 
ci/i^M^ 1 + Ca/^ - ^ - |a;| _Al_2 + o(|ir| _2,i_8 ~ :L ), where ci and C2 are constants. Therefore, the long tails which 
correspond to index [i prevail at large distances and the mean value is infinite. However, in the limit of long time 
the relative contribution to p(x,t) from pi diminishes for a given x since falls faster than fi(t). To demonstrate 
that, we need to estimate the ratio ^1/^2 = /1 ^ / f 2 ^ • That quantity is presented in Fig.l. Its time dependence can 
be very well reproduced by the function 1/t and this pattern is generic for all values of A, /z, and 8. Consequently, 
the contribution from p\, which originates from the Levy process with the order parameter /i, gradually fades away. 
Instead, the term corresponding to fx + 1 dominates the distribution at large time and p\ enters the asymptotic 
expression only with a small weight. 

The transport in superdiffusive systems is used to be characterized by fractional moments of the order 8, (\x\ ), 
where 5 < fi, since all higher moments, in particular the variance, are divergent. That moment is easy to evaluate as 
the Mellin transform \i from the functions pi(x,t): 

poo 

(\x\ s ) = 2 / x 5 p(x,t)dx = 2Nlf^xi(-5-l) + af 2 7 s X 2(-6-l)} 
Jo 

Fig. 2 presents the fractional moments for a few values of the parameter A. The dependence on time is algebraic and 
the power rises with A; for A = the moment is ~ t s /^ +e \ 

Finally, let us mention the case of the attractive potential, A < 0. We require that also a < 0, in order to avoid 
shrinking of p 2 with time (c.f. second equation in Eq. (|16|) L which is unphysical. The negative a, in turn, results 
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FIG. 2: (Colour online). The fractional moments calculated from Eq.fTTJl for S = 0.4, n — 0.5, 6 — 1, and a few values of A. 

in negative p(x, t) for large time, since p\ falls faster with time than p2- Therefore, the above solution of Eq. (ll3p is 
correct for A < only if time is not very large. The fractional moment rises slower than for A = and also weaker 
than algebraically, which is presented in Fig. 2. 

In summary, we have studied the stochastic systems in which the nonhomogeneous structure of the medium is 
reflected not only by an external potential but also directly by the random force in the form of the Levy distribution. 
Such systems are described by the fractional FPE with the variable diffusion coefficient; they have been solved in the 
limit of small wave numbers. The system can have the same stability properties as the driving noise; it is the case for 
the linear drift. However, the other drift we have considered, ~ sgn(x), requires the additional Levy process and then 
the system has the double scaling. That second Levy distribution is characterized by larger order parameter and its 
weight rises with time. 
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